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ABSTRACT 


This paper investigated delamination failure in a unidirectional curved 
composite laminate. The curved laminate failed unstably by delaminations 
developing around the curved region of the laminate at different depths 
through the thickness until virtually all bending stiffness was lost. 
Delamination was assumed to initiate at the location of the highest radial 
stress in the curved region. A closed form curved beam elasticity solution 
and a 2-D finite element analysis (FEA) were conducted to determine this 
location. The variation in the strain energy release rate, G, with 
delamination growth was then determined using the FEA. A strength-based 
failure criterion adequately predicted the interlaminar tension failure which 
caused initial delamination onset. Using the G analysis the delamination was 
predicted to extend into the arm and leg of the laminate, predominantly in 
mode I. As the initial delamination grew around the curved region, the 
maximum radial stress in the newly formed inner sublaminate increased to a 
level sufficient to cause a new delamination to initiate in the sublaminate 
with no increase in applied load. This failure progression was observed 
experimentally . 

INTRODUCTION 

With the increased use of laminated fiber reinforced composite materials 
in primary aircraft structural components, the ability to understand and 
predict their failure modes becomes paramount. Because of the low 
through- the -thickness strengths of composites one of their prime modes of 
failure has been delamination [1]. The high interlaminar stresses which may 
cause delamination can arise from material discontinuities, geometric 


discontinuities, or eccentricities in the load path. One example of a generic 
structural component which has the possibility of all three of these 
interlaminar stress sources is the curved laminate. Some typical uses of a 
curved composite laminate are shown in fig. 1. Mostly, the laminates are of a 
multidirectional layup. The material property mismatch between adjacent plies 
of different orientation may create singular interlaminar stresses at the free 
edge causing edge delamination [2], Another failure mode occurs when the 
membrane stresses in off-axis plies are high and cause transverse tension 
cracking across the width [3]. Mathematically singular interlaminar stresses 
may occur where the matrix crack meets the adjacent plies. These stresses may 
cause delamination at this location. Delamination may also occur as an 
interlaminar tension failure caused by tensile radial stresses created from 
the bending of the curved laminate. 

Several studies in the literature have focused on determining damage in 
curved laminates. In reference 4, a finite element analysis was used to 
determine the strain distribution for the design of curved composite frames 
for the UH60 Black Hawk Helicoptor. In reference 5, a strength-based 
criterion based on the Tsai-Hill criterion was used to predict if an off-axis 
ply failure or an interlaminar tension failure would occur in multidirectional 
curved laminates. In reference 6, a finite element analysis and a closed- form 
solution was conducted to determine the stress distribution in a curved 
composite frame taken from the V22 Osprey. The Tensor Polynomial criterion 
was used to predict failure. In reference 7, an elasticity solution was used 
to determine the stress distribution In a semi-circular composite beam. In 
reference 8, composite angles of different lay-ups were tested and analysed. 
Different failure modes were noted. In one lay-up the failure was 

delamination from an interlaminar tension failure and in another failure was 
matrix cracking in a 90° ply followed by delamination initiating from the 
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matrix crack. A strength-based failure criteria predicted the different 
failure modes but yielded generally conservative predictions. 

A review of some of the work conducted on curved composite beams and 
frames was given in reference 9. Reference 9 notes that predicting 
delamination is difficult to do using strength-based failure criteria because 
adequate test methods for determining the through- the- thickness strengths, 

and r i 2 f> do not exist. A further problem with the use of strength-based 

criteria in the above examples is in the determination of the stresses. 
Although a stress analysis may be useful to identify the areas that are highly 
stressed, a strength-based failure criteria may only be used when no singular 
stresses are present. As soon as some damage or other discontinuity is 
present, such as a radial delamination, a matrix crack, or a free edge, 
stresses cannot be evaluated at the singularities, and a strength criteria 
should not be used. Fracture mechanics based failure criteria offer a 
technique to predict the strength and life of a component with a singularity 
caused by a discontinuity. An attempt to use a fracture mechanics approach to 
predict the static strength of a curved laminate containing both a matrix 
crack and a radial delamination was given in reference 10. The failure loads 
predicted were generally unconservative, possibly because the 2-D finite 
element analysis did not account for the effects of the free edge. In 
reference 11, delamination was assumed to exist in a production part in a 
curved solid laminate segment of a fuel sponson. The failure was assumed to 
occur during manufacture due to residual thermal stresses. A 2-D coupled 
boundary element fracture mechanics approach was used to determine how large 
the delamination could become before the component failed for a given design 
load. It was determined that long before total failure of the component the 
delamination could be easily detected by NDI methods. A fracture mechanics 
based methodology for predicting the static strength and the fatigue life of 
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other structural components which experience delamination initiating from 
singular stress sources has been proposed [12-15]. This methodology uses 
analyses to determine the critical value of strain energy release rate, G , 
and its components along with delamination onset data obtained from generic 
test coupons [16-20] . 

The purpose of this paper is to predict the maximum load a unidirectional 
curved laminate can sustain. Unidirectional laminates were used to isolate 
the interlaminar tension failure, that is it was insured that no edge 
delamination or delamination from cracked plies would occur. Simple 
right-angled coupons were analysed and tested to represent the loading on the 
angle bracket conf iguration shown in fig. 1. First, the location of initial 
failure was determined using a stress analysis assuming an undamaged component 
containing no singular stress sources. A closed form curved beam elasticity 
solution and a finite element analysis (FEA) were used to compute the 
interlaminar stresses. Secondly, a fracture mechanics analysis was conducted 
to determine damage progression after the onset of delamination. This 
characterization was done by determining how G varies with delamination 
growth. A methodology for computing the variation of G for an embedded 
delamination with two delamination fronts was developed. Further damage 
progression was then determined by analyzing the curved laminate with the 
delamination present. 


NOMENCLATURE 
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Width of beam 
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Angle of delamination in clockwise direction 
Displacement in the radial direction 

Displacement in the tangential direction 

Angle around curved region from leg 
Anisotropic parameter for moments 
Length of finite element at delamination front 
Poisson's ratio 
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Radial stress 

Radial stress caused by moments 

Radial stress caused by end loads 
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MATERIALS 


Twenty- four ply, graphite/epoxy (AS4/3501-6) unidirectional curved 
laminates were manufactured at NASA Langley Research Center. The laminates 
were cured according to the manufacturer's instructions. The curved laminates 
were layed up over a solid aluminum tool in panels 300 mm wide to give a 90° 
bend. Following curing, the panels were cut into 25 mm wide specimens. The 
final dimensions of the curved laminates are given in figure 2. Double 
Cantilever Beam, DCB, specimens were also manufactured. The DCB specimens 
were 24-ply unidirectional coupons with a sheet of 0.127 mm (0.5mil) Kapton 
film placed between the 12th and 13th plies prior to curing to simulate a 
delaraination. The dimensions of the DCB specimens were given in reference 17. 
The following values of fracture toughness were obtained from delamination 
growth initiating at the insert: 

69.0 < G Ic < 99.0 J/m 2 (mean - 80.0 J/m 2 ) 

Values for transverse tensile strength obtained from flat specimens were taken 
from reference 21 and were: 

36 . 5 < a ^ i 69.0 MPa (mean - 53.8 MPa) 

Specimens to determine the elastic properties of the composite were also 
manufactured. The elastic tensile and shear properties used in the FEA and 
the closed form elasticity solution were obtained using ASTM standards D3039 
and D3518, respectively. The following average properties from at least three 
tests were obtained assuming transverse isotropy: 

E u - E. - 138 GPa E 22" E 33“ E r“ 10 . 0 GPa 

G 12“ G rr 4l47 GPa "12 “"13" °‘ 3 "2l“"3l“"rr 0 0218 
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The average volume fraction, V^, of the specimens was determined from 
FAWxNxlOO 


- FDxt — where FAW is the fiber areal weight, N is the number of plies, 

FD is the fiber density, and t is the thickness of the specimen. For the DCB 

specimen, V^. was determined to be 55.9 percent. In the curved laminate 

specimens, in the curved region was 54.7 percent and in the leg was 56.1 

percent. In the specimens used to obtain elastic moduli, - 59.1 percent. 


EXPERIMENTAL PROCEDURE AND RESULTS 


A special fixture for testing the curved laminates in a standard tensile 
testing machine was manufactured and is shown in fig. 3. A 890 N (200 lb) 
load cell was incorporated into the test fixture to monitor accurately the 
small loads associated with matrix failures in composites. The load was 
applied to the specimen at a distance L of 50 mm from the curved region/arm 
intersection via a hinge clamped to the specimen. The other end of the hinge 
was held in a clamp that was free to pivot. The displacements were small 
during the test because of the high stiffness of the unidirectional curved 
laminate. Therefore, the load applied via the hinge was considered to remain 
vertical during testing. A displacement transducer, not shown in fig. 3, was 
attached to the side of the loading arrangement to measure vertical 
displacements during the tests. The edges of the specimens were painted with 
a brittle white paint to aid visualization of interlaminar failures. The 
load and displacement were monitored on an X-Y plotter to determine the 
failure load. The laminates were loaded quasi- statically in load control. 

The maximum static loads per unit width for damage initiation in the 
unidirectional curved laminates are shown In table 1. A significant scatter 
can be seen in the results. A similar scatter was reported in reference 6 and 
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is consistent with the large scatter in transverse tensile strength in 
reference 21. The specimens were labeled sequentially from the location they 
were cut from the panels. There was no trend found between the strength and 
the location where the curved laminate was cut from the panel. The normalized 
radius location of the first delamination is also given in table 1. The 
damage in three different specimens is shown in fig. 5. Damage initiated, 
fig. 5a, as a single delamination around the curved region. Growth was 

initially unstable and extended into the arm and the leg. Two delaminations 
around the curved region are shown in fig. 5b. the lower delamination is at an 
(R-R. )/t of 0.23. Final damage was an accumulation of delaminations at 

different locations in the thickness of the curved region until the specimen 

lost virtually all bending stiffness. The longest delamination in fig. 5c is 

at an (R-R. )/t of 0.46 and was presumably the first delamination to occur, 

inner ' 

For the upper two photographs shown in fig. 5 to be taken, the load was 
reduced manually when audible accoustic emissions were heard. 


ANALYSIS 


The following section details the analysis required to characterize the 
onset and growth of a delamination in the unidirectional curved laminate. The 
location of the highest radial stress, a to cause delamination was determined 

using a closed form elasticity solution for an anisotropic curved beam. A 2-D 
FEA was conducted to verify the beam solution and determine the variation of 
G with delamination growth from the location of highest radial stress. 

Closed Forra Anisotropic Elasticity Solution for a Curved Beam 

The loading on the curved laminate used in the experimental work may be 
reduced to a curved beam with an end load, P, and moment, H, at the free end 
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of a quarter of a curved circular beam, figure 4. The stresses created In the 
beam may be calculated by summing the stresses caused by the end load and the 
moment. A closed form elasticity solution for the stresses in a curved beam 
with cylindrical anisotropy and subjected to end forces and moments was given 
in reference 22 and is repeated in Appendix A using the notation give in fig. 
6. It was assumed that the location of the maximum radial stress in the 
curved portion may indicate where delamination onset will occur. The maximum 
radial stress was determined by solving 



dr 



de 


p 

for R and 9 , where o is the radial stress created by the end load P and 
' r 

Is the radial stress created by the moment M. The solution yielded a 
r 

maximum value of a ^ 4 - at R- 6.21 mm and 0=“O° . This value of R was 

r r 

approximated to R=6.25 mm which was the closest radius to two ply Interfaces. 
However, the boundary conditions are for a curved beam and not a laminate with 
straight attachments. To determine the effect of straight legs attached to a 
quarter circle a finite element analysis, described in the next section, was 
conducted . 

Finite Element Analysis 

A 2-D FEA was conducted using the finite element code NASTRAN [23] using 
two different models. Model 1, shown in fig. 7a, was used to determine the 
stresses in the curved region of the beam assuming no damage in the part. The 
model had one element per ply thickness and one element per one degree sweep 
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around the curved region. The model had a unit end load per width (1.0 N/mm) 
applied at a length of 50 mm along the arm to simulate the experimental 
conditions. Model 2 with a refined mesh, fig. 7b, was used so that the 
virtual crack closure technique, VCCT, [24] could be employed to determine 
total G and its components corresponding to delamination growth. The mesh was 
refined at R*=6.25 mm which was the radius determined from the closed form 
elasticity solution (and the FEA as described later) where delamination was 
expected to initiate. Near a radius of R-6.25 mm the finite elements measured 
h/4 by h/4.6 where h is one ply thickness. At a radius of R=6.25 mm a free 
surface was included in the model by the use of coincident nodes to represent 
the delamination. The free surface extended around the complete 90° of the 
curved region. The coincident nodes were restrained together using 
multi-point constraints (MPCs) . By releasing different MPCs in different 
analysis cases, several delamination lengths could be modeled during one FEA 
run yielding an efficient computational technique. To reduce the number of 
elements in the model the length of the arm was shortened to L^=5 mm. A load 

and moment were applied to the shortened arm, fig. 4. The load and moment 
were equivalent to an end load of 1 N/mm applied at a moment arm length of 
50 mm. 

Calculation of Strain Energy Release Rate 

The virtual crack closure technique was used to obtain total and 
individual modal values of G at the delamination fronts. The individual 
modes of G were calculated at the delamination front using the following 
equations in polar coordinate notation, fig, 8, 

gY t - 1 F 4 f 5 B - s c ] m 
1 TTZ r l r r J 
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1 


( 4 ) 



2 A w 


F* 



) 


where and are the radial and tangential forces at node A respectively, 

A is the length of the finite element at the delamination front, 5 ^ and 6 $ 

are the radial and tangential displacements at the nodes indicated by the 

superscripts, and w is the width of the model. Total G calculated by VCCT, 

VT . VT VT 

G , was simply a summation of G^ and 


An alternative method of calculating total G to check the VCCT results 
was also used. This method involved determining the difference in total 
strain energy from two finite element runs at two delamination lengths, 1^ and 


! 1+ 1 » thus 


g se - 


2 w (1 - l i+1 > 


S (P «). - S (P 6 ) 

l i+l 


(5) 


where 2(P 6 ) is the sum of the product of the applied loads and the 
displacements (in the loading direction) of the loaded nodes. The resulting 


C IT 

G approximately applies at a delamination length of (l^+l^ + ^)/2, 


ANALYTICAL RESULTS 

Stress Distribution 

Figures 9 and 10 show the radial, transverse shear and tangential 
stresses, respectively, as calculated by the closed form elasticity solution 
and FEA. The stresses were normalized by P/w. The stresses are plotted 
versus the normalized radius at an angle 9 — 25° . The choice of this angle 
will become apparent in the next paragraph. In fig. 9 the agreement for the 
radial and transverse shear stress between the two analytical methods was good 
and confirms that the radial stress was highest at (R-R^ nner )/t “ 0.42 
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(R-6.25 nun). The transverse shear stress was approximately an order of 
magnitude smaller than the radial stress. The agreement in fig. 10 was also 
good and shows that the laminate has a maximum tensile tangential stress at 
the inside edge of the curved region. 

Figure 11 shows the normalized radial stress plotted around the radius at 
R-6.25 mm. The closed form elasticity solution agreed well with the 

FEA for 2O°<0<6O° . However, as 0 approaches the boundaries the two solutions 
diverge as expected from St. Venant's principle. According to the FEA the 
radial stress reaches a maximum at 0~25°. Therefore, the location for the 
onset of delamination was determined from figures 9-11 to be at R=6.25 mm at 
an angle 0 — 25°, 

Strain Energy Release Rate Variation with Delanination Growth 

Next, the variation of G with delamination growth was determined. The 
delamination had two fronts which may have different values of G. Referring 
to fig. 6, the delamination front may grow counter-clockwise or clockwise from 

the initiation point (R*=6.25 mm, 0=25°) as indicated by the letter a or b, 

2 r q 

respectively. The delamination length is given by a+b where a - — — and 


b 


2 7T R B 
360 


and a and are measured in degrees . 


2 2 

Figure 12 shows the normalized G distribution (G w /P ) with delamination 
growth in the a direction only, that is holding 0 equal to zero. The 
delamination was predominantly mode I around the curved region. The 


. VT VT VT 

normalized values of G obtained by VCCT (G^. + G^) closely agree with 

SE 

values obtained using the change in global strain energy technique, G . Both 
SE VT 

G and G^ increase to a maximum at a=22° and then reduce. At a-25 0 the 


delamination is at the curved region/leg boundary. 
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Figure 13 shows the normalized G distribution in the f) direction, holding 
a equal to zero. Again the delamination growth was predominantly mode I 
around the curved region. The total G results from VCCT and the change in 

VT SE 

global strain energy technique again agree well. Both G^ and G increase 

VT 

from zero to a maximum and then decrease. For G^ the maximum is at /?=48° 
SE 

while for G the maximum is at /3-50 0 . The difference between the two peaks 

VT VT 

is caused by the increase in Fig. 13 shows that G ^ increases 

throughout the curved portion of the beam. However, as the delamination grows 
into the arm, the horizontal forces caused by the moment would decrease. 
VT 

Hence, G ^ should reach a maximum at the start of the arm, and decrease as 
the delamination grows into the arm. 

Figures 12 and 13 show how G varies when either one of the two fronts 
grows individually. However, they do not show how the two fronts grow 

simultaneously. Consider an angle <f > , where <£-**+£. Also, assume that the 

delamination grows equally in the a and b direction, hence, a-/?- ~ 2 ~* The 

result of assuming this delamination growth behavior is given in fig. 14. 
Note that when ^=50°, the delamination has grown fully around the curved 
region to the leg but has not yet reached the arm, fig. 6. For small angles 
of <j>, G a and G^ are similar and the delamination should grow equally in the 

a and fi direction. However, at larger angles of <f> , G becomes less than G^ 

because is decreasing more rapidly in the a direction than in the 

direction, fig. 11. Hence, the delamination would continue to grow in the £ 
direction and not in the a direction. Thus, any subsequent data points in 
fig. 14 are invalid. Delamination growth must be represented incrementally 
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and and G^ compared at each increment. If G at one end of the 

delamination is higher than that at the other, then the delamination will grow 
in that direction only. Following that incremental growth another analysis 
must be conducted, and so on until the delamination arrests or failure occurs. 
Hence, a methodology was developed to ensure a constant G existed at both 
delamination fronts. A criterion was employed that if G was two percent 
greater at one delamination front than at the other, then growth would occur 
in that direction only. If the difference was less than two percent then the 
delamination would grow in both directions simultaneously. For the first 
increment a 5° step was taken from the Initiation site. From there, 
incremental delamination growth steps of 2.5° or 0.273 mm were assumed. The 
results are shown in fig. 15. Note that G is plotted against a which ranges 

from 0 to 25° whereas G^ is plotted against which ranges from 0 to 65°. 

From the onset location, the delamination grows simultaneously in both 
directions as Indicated in figure 14. However, at the third increment G^ was 

3.7 percent higher than G . Hence, delamination growth was assumed to occur 

a 

in the direction alone, and is shown in fig. 15 by the next G value being 

a 

plotted directly above the previous one at the same value of a. At this 
increment G^ and G^ have less than a 2 percent difference. Hence, 

delamination growth will occur in both directions simultaneously shown by the 
next G^ and G^ values being plotted at an increased value of a and 

respectively. This procedure was continued until the delamination reached the 
leg. 

SE 

A plot of normalized G versus <f> is shown in fig. 16 using the results 
of fig. 15. This plot details the variation in normalized G with 

delamination growth taking into account the two delamination fronts. The 


14 



delamination grows unstably until it has grown 60° around the curve. Total G 
then decreases and the growth becomes stable. Fig. 17 shows the variation of 


the Gj/G ratio 

as 

the 

delamination 

grows . Delamination 

growth is 

predominantly mode 

I. 

with 

pure mode I 

delamination growth 

at small 


delamination lengths. 

Damage Progression in the Unidirectional Curved Laminates 

Fig. 18 shows the normalized radial stresses plotted versus the 
normalized radius of the curved laminate with a delamination present at 
(R-R^ )/t-0 . 42 , Three specific delamination lengths taken from fig. 16 are 

shown. As the delamination grows around the curved region the maximum radial 

stress in the newly forming inner and outer sub-laminates increases. At <£=90° 
G is decreasing (fig. 16) and the maximum radial stress in the inner 
sub- laminate has reached a value similar to that responsible for the formation 
of the original delamination at R=6.25 mm as shown in fig. 9. Therefore, a 
new delamination should start in the inner sub- laminate as was observed 
experimentally, fig. 5. No further analysis of delamination growth from this 
location was conducted but it was assumed that this phenomenon would repeat 
sequentially and result in the final failure in fig. 5. 

FAILURE LOAD PREDICTIONS AND DISCUSSION 

To determine the maximum loads required to initiate damage in the 
unidirectional curved laminates, a suitable failure criterion must be used. 
Because of the high tensile strength in the fiber direction it is not 
necessary to consider the a $ stresses. Also, the values of are an order 

of magnitude smaller than the a stresses, and may be considered negligible. 

Therefore, it is sufficient for a strength approximation to assume that 
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delamination occurs when a where is the interlaminar normal 

r 3f 3f 

strength. However, as noted in the introduction, no suitable test for 
determining interlaminar normal strength has been developed. Hence, 
transverse isotropy must be considered and a ^ approximated by From 

Fig. 9 the maximum normalized radial stress was given as (a^ w /P) - 4.37 mm’ 1 
occurring at an (^"^^ nner )/ t =0.42. Therefore, a range of failure loads per 
unit width can be determined: 


w 


u 2f 

i° x w/P) 


8.36 < <15.8 N/mm (mean - 12.3 N/mm) 


w 


Table 1 shows that the agreement between the range of failure loads and the 
range of predicted loads was good. The mean failure load was over predicted 
by 15 percent. The predicted location of the delamination through the 
thickness was similar to that observed in the experimental work, Table 1. 

The plot of normalized G versus <j> shown in fig. 16 shows G beginning to 
decrease at <^60° , This type of decrease in G with delamination growth is 
often interpreted as an indication of stable delamination growth. However, 
the delamination usually grew well into the arm and leg, fig. 5. Delamination 
growth was shown in fig. 17 to be largely mode I dominated. Therefore, using 
the criteria that G-G c «Gj c at failure a predicted load per unit width to 

further extend the delamination may be calculated. From fig. 16 at <£~90°, 


(G w 2 /P 2 ) -1.60 E-2 mm/N: 




2.01 < < 2.41 N/mm (mean - 2.24 N/mm) 

w 


where G^ was obtained from the DCB specimens. These values of applied load 
per unit width are well below the load already applied to the curved laminate, 
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indicating that the delamination will extend into the arm and the leg until G 

decreases below G . 

c 

At a delamination length of ^=90° the maximum normalized radial stress in 

the inner sublaminate from fig. 18 was (a w/P)- 4.10 N/mm ^ at an 

(R-R. )/t ~ 0.20. A new delamination will occur at a load per unit width 

inner 

between 8.90 and 16.8 N/mm. This range of loads is approximately equal to the 
range of loads already on the specimen. Also, the predicted location of the 
second delamination is similar to that shown in fig. 5b. Hence, damage 

accumulating at different locations through the thickness would be expected to 
be unstable, as observed experimentally. 


SUMMARY 

This paper investigated delamination failure in curved composite 
laminates, A unidirectional curved laminate was tested quasi -statically and 
analyzed. These laminates failed unstably by delaminations radiating around 
the curved region of the laminate and accumulating at different locations 
through the thickness until most of the bending stiffness was lost. 
Delamination was assumed to initiate at the location of the highest radial 
stress within the curved region. A closed form curved beam elasticity 
solution and a finite element analysis (FEA) were used to determine the 
location of the highest radial stress. A strength-based failure criteria 
predicted the failure load for the onset of the first delamination. The FEA 
was then used to determine the change in strain energy release rate, G, with 
delamination growth in both directions from the initiation site. The strain 
energy release rate was found to be predominantly mode I. Strain energy 
release rate initially increased with delamination growth and then decreased 
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as the delamination reached the straight arm and leg. However, the value of 
G with delamination fully extended around the curved region was larger than 
the critical value of G determined from Double Cantilever Beam specimens. 
Thus, delamination growth was predicted to extend into the arm and leg of the 
laminate as observed experimentally. As the initial delamination grew around 
the curved region, the maximum radial stress in the newly forming inner 
sublaminate increased to a sufficient level to cause a new delamination to 
initiate in the sublaminate with no increase in applied load. This failure 
sequence was also observed experimentally. 
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Appendix A - Closed Form Elasticity Solution for the Stresses in an 

Anisotropic Curved Beam 


Reference 22 gave the elasticity solutions for the radial, tangential, 
and shear stresses in a curved beam subjected to pure bending and an end 
force. The beam was assumed to possess cylindrical anisotropy and its pole is 
located at the center of the arcs which form the beam contour. The 
expressions are repeated here with the notation given in the Nomenclature and 
in fig. 6 
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Panel No. - Specimen No. 


Maximum load 


(R -R. ) 

inner 



per unit width 
(N/mm) 

of first delamination 

6-4 

7.62 

0.49 

6-8 

14.60 

0.48 

7-1 

8.68 

0.46 

7-2 

10.05 

0.48 

7-6 

13.36 

0.48 

7-7 

13.09 

0.35 

7-8 

7.23 

0.33 


f 

7.23 < ~ < 14.60 N/mm 

Mean = 10.66 N/mm 


Table 1 Experimental Failure Loads in the Curved Laminates 
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Curved laminate configurations and failures. 
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Dimensions of curved laminate. 
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- Curved laminate test fixture. 




Experimental Elasticity solution Finite element 
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Damage in unidirectional curved laminate. 


Location of 
delamination 
s onset — \ 



30 


Fig. 6. - Notation for unidirectional clurved laminate. 
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Virtual crack closure notation. 
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Interlaminar stresses through thickness. 
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Fig. 10. - Tangential stress through thickness. 




11. - Radial stress around curved region. 
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Fig. 14. - G Variation with delamination growth in the a and 3 
direction. 
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Fig. 17. - G|/G Ratio with delamination growth. 
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Fig. 18. - Radial stress through thickness in presence of a 
delamination. 
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